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Abstract 

A vrvr, KK, and pp{ujuj) fully coupled channel model is used to predict the 



scattering from threshold to 2.0 GeV. As input the S-matrix is required to 
exhibit poles corresponding to the meson resonance table of the Particle Data 
Group. As expected, the vrvr inelasticity is very strongly related to the open- 



^ . ing of the KK channel near 1 GeV, and the opening of pp{Att) and ujlo{Qtt) 

channels in the 1.5 GeV region. The predictions of this model are compared 
to the various elastic vrvr tttt amplitudes, that were obtained from analyses 
of vr^ p vr^vr^ n data. The role of the various resonances, in particular the 
glueball candidate /o(1500) and the /j(1710) is investigated. 
PACS:13.75.Lb; 14.40.-n; 13.25.-k; 12.39.Pn 
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I. INTRODUCTION 



In heavy ion collisions and in anti-nucleon physics many pions are produced in the final 
state. Pion-pion scattering therefore plays an important role in the final state interaction of 
these processes. Our knowledge of tttt scattering is incomplete, in particular above M^^^^ ^ 
1 GeV. The dynamics of vrvr scattering is often described by effective meson-exchange in the 
t-channel. The mechanism of t-channel exchange works very well in, for example, nucleon- 
nucleon scattering. However, for vrvr scattering one may also consider the presence of the 
many resonances in the s-channel, a feature that is typical for tttt scattering but not for 
nucleon-nucleon scattering. A list of relevant meson resonances with their properties can be 
obtained from the compilation of the Particle Data Group For energies larger than 1 
GeV, an additional aspect is the strong coupling of the vrvr channel to KK and multi-pion 
channels, for example, four and six pions. 

Previous models for S-wave tttt scattering have been applied from threshold to 1.4 GeV. 
For example, separable potential models have been considered in an approach in which 
the major focus was to obtain the proper low energy behavior of the tttt S-wave phase shift, 
using the right combination of attraction and repulsion in the diagonal vrvr interaction, and 
also to better understand the structure of the meson resonance /o(980). More sophisticated 
one-meson exchange models ||^-^ have been used up to 2 GeV, still within the framework 
of coupling vrvr and KK channels. 

In an earlier paper we reported on results within a coupled three-channel model for 
S-wave scattering in the 1.0 - 2.0 GeV energy region, and we were encouraged by its results 
when compared to the then available S-wave data. The relevant experimental vrvr phase 
shifts and inelasticities, specifically vr"'"vr~, can be found from the analysis by Protopopescu 
et al [|10| , Grayer et al |jlll , Hyams et al fl^ , Bugg et al |jl3| and Kaminsky et al |]14| . These 
groups have extracted the pion-pion scattering amplitudes in the M^t^ = 0.60 - 1.78 GeV 
region by obtaining them from an analysis of the reaction vr~ p — * vr^vr^ n. Kaminsky et 
al included also data from the same reaction with a polarized proton target, which allowed 
them to separate the contributions due to vr and ai exchange. The low energy vrvr data were 
published by Ref. |TB|. From all these extracted data for the reaction vr+vr^ vr+vr^ it is 



apparent that there is a rather strong energy dependence of the phase shifts as well as the 
inelasticities for S, P, D and F-wave vr"''vr~ elastic scattering. We do not know of any data 
above 1.78 GeV, so a model reproducing the existing data may eventually also serve as a 
basis for extrapolation to higher energies. Such an extrapolation would be required to study 
for example the final state interaction in the process pp —>■ vr+vr^ [|T6|. 

In view of the presence of an extensive set of resonances as given by the Particle Data 
Group we have opted in this paper for inclusion of the various known resonances in 
our model for vrvr scattering. One can then address the question, what is the role of these 
resonances for elastic vrvr scattering in the various partial waves in the region of 0.3 - 2.0 
GeV. In the energy range up to 2.0 GeV the dominant channels are vrvr, KK, pp and uu. 
Within our model, as no spin effects have been taken into account, the pp and ujuj channels 
can be described by a single effective channel, e.g. pp. 

For example, the set of S-wave resonances that play a role for the vrvr channel are 
the /o(980) and /o(1370) resonances, but one should also include the recently discovered 
/o(1500). While many conjectures about the nature of the /o(1500) have been made, its 
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role for tttt scattering is still an open question, and therefore will be one of the questions 
addressed in this paper. Another question is the J value of the reported /j(1710) reso- 
nance. So far one only agrees that J should be even [|I|. In this approach one can study 
the implications of the presence of an isoscalar /o(1710) resonance for S-wave vrvr scattering, 
as well as the implications of the existence of an isoscalar /2(1710) resonance for D-wave 
TTTT scattering. In practice we will impose the condition that the scattering matrix, corre- 
sponding to a particular angular momentum, should have poles at the complex energies of 
known resonances with appropriate quantum numbers and, from this input, one can study 
the implications for tttt scattering. 

The model that we will describe, can be applied for any angular momentum. In this paper 
it will be applied for S, P, D and F- waves. For each angular momentum the important input 
will be the location of the corresponding resonances as given by the Particle Data Group [|1| . 
Sections II and III describe the model and show the analytic expressions of the S-matrix for 
the various angular momenta. In sections IV and V the results for respectively / = and / 
= 2 are discussed. This order of angular momenta is chosen because the resonance /j(1710) 
may either contribute to the S-wave or D-wave scattering. Then the results for / = 1 and / 
= 3 are described in sections VI and VII respectively. Finally a discussion of all results and 
conclusions follows in section VIII. 



II. MODEL 

As in Ref. we consider three coupled channels of tttt, KK and pp, and label them 
respectively with the channel index z = 1, 2, 3. Our method of derivation of the scatter- 
ing amplitude is generalized to all angular momenta. For each angular momentum / the 
Lippmann-Schwinger equation for the zj-element of the T-matrix at total invariant energy 
squared s, is 

<p\Tl^{s)\q >=<p\Vl^\q > "^E / dp'p" < pK\p' > G,{p',s) < p'\Tl.{s)\q > . (1) 

Here Gk{p' , s) is the propagator for channel k and V^^ is part of the potential. In this case 
the potential as well as the T-matrix are both 3x3 matrices. 

We choose a separable form of the potential matrix elements for each angular momen- 
tum / 

<p\Vy\q>=g\{p)\\^g]{q). (2) 
As a consequence the T-matrix elements for each / are also separable 

<p\Tl^{s)\q>= g\{p)Tl^{s)g]{q), (3) 

where rj, satisfies 

''J 

rUs) = X\,-j:X\,Xl{s)Ti^{s). (4) 
fc=i 

The functions X[ are defined in terms of the propagator Gk and vertex functions gl 
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Xlis) = ^ / dp'p"gl{p')G,{p', s)gWl (5) 

and the functions T?j(s) can be obtained in closed form. 

In the following we will suppress as much as possible the label I of the angular momentum. 
For example if the matrix A for a given angular momentum / is defined by 

-^ij = \j + (6) 
the element th has the general form 

(1 + A22-^2)(l + A33X3) — A23X3A32X2 , . 

= d^) • 

The elastic vrvr scattering amplitude for angular momentum / is then given by the corre- 
sponding element Tij with i = 1 and j = 1, of the T-matrix of Eq. (^. The model has a 
resonance of angular momentum I if the corresponding T-matrix has a pole in the complex 
energy plane at an energy whose real part is the resonance mass and whose imaginary part is 
half the resonance total width. Such a pole corresponds to a zero of det{A) at that complex 
energy, as can be seen from Eq. (|^). The full expression for det{A) is 

detiA) = (1 + AnXi)(l + A22X2)(1 + A33X3) - A23X3A32X2(1 + XnXi) 

— Ai2X2A2lXi(l + A33X3) — Ai3X3A3iXi(l + A22X2) 

+ (Al3A32A2l + Ai2A23A3i)XiX2X3, (8) 

with the couplings A'^- satisfying 

= (9) 
We choose the form factor gl{p) of Eq. (j^) to be 

gl((3,p) = .^ f (10) 

where (3 carries an index i as well as /. For the propagator Gi{p', s) we take the form 

G.{p',s)= ^; . , (11) 

P — Pi — 

where s is related to pj by 

s = Aip^ + m^). (12) 

The functions for / = to 3, which are dependent on the range parameters Pi, become 
in this case 
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(14) 



Again the indices i and / of (3 and ^ of p are suppressed. The above expressions follow from 
evaluation of the integral in Eq. for real values of p using the definitions Eq. (|1^) and 
Eq. (0). Details of the evaluation are given in part A of the Appendix. The integral in 
Eq. d^) is defined for Re(p) > 0, and is dis-continuous over the unitarity cut for positive real 
p. Therefore it cannot be analytically continued for complex p with Im(^) < 0. However, 
the expressions of Eqs. ([T^), (0), (|T^), (0) are continuous and analytic for all complex 
momenta p, except for a pole on the negative imaginary axis, and can therefore be used 
as analytic continuation of X\{f3,p,s) to complex values of p, even to the resonance region 
where Im(]?) < 0. One also notes the familiar property of Xl{P,p, s) 

Xl{(3,p,s) = [Xl{(3,-p\s*)r. (17) 

As mentioned previously, the channels 2 = 1, 2, 3 correspond to vrvr, KK and pp respec- 
tively, and we use mi = 0.1396 GeV, m2 = 0.4937 GeV and ms = 0.7680 GeV. Since no 
spin aspects are considered and the mass of the u meson is very close to the mass of the 
p meson, the pp channel can effectively be viewed as representing also the uu channel for 
angular momenta where 1 = 0. We ignore the width of the p meson, and assume that Air 
and Gtt channels are dominated by pp and uju respectively. The relative momenta, pi, in the 
i*'* channel, are related by the kinematic condition 

I = = pI + ml = pI + ml = pI + ml (18) 

This equation defines the invariant mass M^^^^. For / = the form factor in coordinate space 
has the familiar form 

9^^^) = (19) 

The above method is equivalent to previously given expressions for the S-wave elastic 
TTTT scattering amplitude, 

_ D{-pi,p2,P3) . . 

D{pi,P2,P3) 

where Su is the vrvr S-matrix element. D is the Jost function, and D = det{A), i.e. 

D{pi,P2,P3) = [{Rl + Al)(i?2 + A2)(i?3 + A3) - AURs + A3) - A?3(i?2 + A2) 

-A-URi + Ai) + 2Ai2Ai3A23]/(i?i/22i?3). (21) 



5 



The couplings A' and K\j as well as the functions R\ are dimensionless and their relation to 
\\j and Xl of Eq. (j^) are given in Table |. Since the function X' has the property given in 
Eq. (p!?!), the Jost function satisfies 

^(Pi,P2,P3) = [D{-pl, -Pi -Pi)]*. (22) 

For each angular momentum I = J we now require the S-matrix to have poles at the 
known J'^ resonances. A summary of the considered resonances appears in Table ||. For 
example, for J'^ = 0+ the resonances are /o(980), /o(1370) and /o(1500). For each angular 
momentum / = J there are three range parameters and six couplings Xij. The three poles 
are then chosen at the complex energies, {R.eM.^.,r, ImM-^T^) corresponding respectively to the 
mass of the resonance and half its total decay width. Again for = 0^, these energies 
would be respectively (0.980, -0.0250), (1.300, -0.200) and (1.503, -0.060) GeV. We require 
det{A) of Eq. (|]) to have zeros at these complex energies, which leads to six constraints on 
the nine parameters. There are several parameter sets that satisfy the above constraints, as 
will be discussed below. 

III. ADDITIONAL FORMALISM FOR J'^ = 0+, / = 

The case J'^ = 0+ needs special attention. We have before pubhshed results ^ for a three 
coupled channel model as described in Section II. One obtains impressive results for phase 
shifts as well as inelasticities in the energy region 0.9 - 2.0 GeV. However, within the model 
described in Section II it is not possible to obtain a realistic description at energies below 
0.9 GeV. The reason is that the potential for the vrvr channel for the S-channel requires 
an attractive term as well as a repulsive one. It has been earlier demonstrated ||^ that 
addition of a second term in the S-wave tttt potential is successful in describing the low 
energy behavior of the / = scattering. We develop therefore a formalism, described in 
Appendix B, adding such an attractive term. In this paper this formalism is only applied 
for / = scattering. 

IV. RESULTS FOR = 0+, 1 = 

Since our previous publication of / = results 0, the Particle Data Group [0] has 
significantly revised its compilation of J^{I) = 0^(0) resonances. At present it contains the 
resonances /o(980), /o(1370), /o(1500) and possibly /o(1710). They occur respectively at 
the complex energies (0.980, -0.025), (1.300, -0.200), (1.503,-0.060) and (1.697, -0.088) GeV. 
However, the real as well as the imaginary part of the pole positions have uncertainties and 
those are rather large for the /o(1370) resonance. 

We are particularly interested in the role played by the "new" resonance /o(1500), which 
has been considered as a glueball candidate. Therefore we start by taking the first three res- 
onances, /o(980), /o(1370) and /o(1500), as input for our model and study the consequences 
for the resulting tttt scattering parameters. These three J = resonances together will im- 
pose constraints on the analytic expression for the S-matrix for / = 0. The S-matrix must 
have three poles at the corresponding three complex energies. These constraints determine 
or restrict the possible values of the model parameters /3o, /32, /^a, Aq, Ai, A2, A3, A12, A13 
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and A23, where, as stated in Appendix B, the index labels the second potential term of 
the TTTT channel. Using these parameter values, one can then immediately determine the 
S-submatrix Sn, describing tttt — > tttt scattering. Since 5*11 for angular momentum / is 
parametrized by a phase shift 6 and an inelasticity rj, viz. 

Su = r/e2^^ (23) 

one obtains the values of 6 and rj as a function of energy. The / = scattering parameters 
are then compared to experimental data. 

There is still an additional freedom in our model. This is related to the location of 
the resonance in the complex energy plane and to the analytic structure of the S-matrix 
for the three coupled channels tttt, KK and pp. Since the relation between energy and 
momentum is quadratic, for each complex resonance energy there are two corresponding 
complex momenta Pi for each channel, one with Im pi < and a second with Im pi > 0. 
This defines two possible sheets in the complex energy plane for each channel. Because there 
are three channels, there are all together eight sheets. The exact location off all resonance 
poles with respect to these eight sheets determines the specific solution for the S-matrix. 

In our standard choice of complex channel momenta Pi for each resonance, the sign of 
Im Pi is such that the momenta Pi are all close to the physical scattering region. The physical 
scattering region has positive real momenta pi for the open channels, and positive imaginary 
momenta pi for the closed channels. This means that a standard resonance corresponds to 
a channel momentum pi with Im pj < if channel i is open and Im pj > if channel i is 
closed. The tttt channel is always open, the KK channel opens at 2mK and the pp channel 
opens at 2 rrip. 

If we put all three resonances /o(980), /o(1370) and /o(1500), on the standard tttt , KK 
and pp sheets, i.e. near the physical region, (where all three channel momenta pi,P2,P3 
have the required imaginary part, in agreement with each channel being open or closed), a 
solution is found for the parameter set (/3j, Aj, Aij) that is given in the first line of Table. p| . 
However, the corresponding prediction for the tttt phase shifts is considerably lower than 
the data suggests, as can be seen by the dashed curve in Fig. Therefore we next consider 
cases where one of the resonances is not close to the physical scattering region. 

From the decay properties of the first three J = resonances one knows that /o(980) and 
/o(1370) both have significant couplings to the tttt and KK channels. On the other hand 
the /o(1500) resonance has a preference to decay into rji]', 7]7], 4tt°, 2tt° and 27r+27r^. It is 
significant that the decay of /o(1500) into tt+tt^ and KK may be rather small. It suggests 
that the /o(1500) resonance could be located on the tt~^tt~ or KK non-standard sheets (i.e. 
its pole location may correspond to a channel momentum pi with a positive imaginary part). 

If we choose to put /o(1500) on the non-standard vr+vr" sheet, (i.e. treat it as a virtual 
resonance, where the relative vrvr momentum has a positive imaginary part), there is a 
significant improvement of the model prediction when compared to the experimental data. 
The predictions for the tttt phase shifts and inelasticities are shown as the solid curve in 
Fig. and Fig. |^. The corresponding parameter set is given in the second line of Table |T|. 
The behavior of phase shift and inelasticity above 1.6 GeV will be influenced further by 
higher resonances. On the other hand, if we put the /o(1500) resonance on the sheet, where 
also the relative KK momentum has a positive imaginary part, the vrvr phase shift obtains 
a significant structure in the 1.2 - 1.4 GeV region. Such a structure is not present in the 
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experimental data. We therefore conclude that the data for vr"*"??" scattering are consistent 
with a /o(1500) resonance that is very different from the /o(980) and /o(1370) and that it 
couples only weakly to tt+tt". At the same time these data do not require its coupling to 
KK to be also weak. 

The next question is if these data are consistent with the existence of a /o(1710) res- 
onance. Indeed it is possible with resonance poles at the complex masses (1.200, -0.250) 
and (1.697, -0.300) to obtain an excellent prediction for the phase shift in the 1.0 - 1.5 GeV 
region that is in good agreement with the experimental data. In Fig. ^ and Fig. ^ the S-wave 
phase shifts and S-wave inelasticities for three parameter sets, are shown when /o(1710) is 
included. In this case one requires a large width for /o(1710) as well as a low mass for the 
/o(1370). The corresponding three sets of parameters are given in the third to fifth line of 
Table The phase shifts in this case (see Fig. ^ show an energy dependence that is closer 
to the data than the model prediction in Fig. |I[ This makes a strong case in favor of J = 
for /j(1710). The possibility of the existence of /2(1710) still remains to be considered and 
the results for J = 2 will be discussed in the next section. 



V. RESULTS FOR J'^ = 2+, / = 

The Particle Data Group lists several J'^{I) = 2+(0) resonances, /2(1270), /2(1525) and 
/2(2010) with the possibility of an additional /2(1710). The first three resonances occur 
at complex energies (1.275, -0.092), (1.525, -0.038) and (2.011, -0.101) GeV. However, a 
resonance with even J has also been observed at (1.697, -0.088) GeV and we will therefore 
follow two avenues. First, in case A we impose on our model the existence of three resonances 
in J = 2 at complex energies (1.275, -0.092), (1.525, -0.038) and (1.607, -0.088) GeV, while 
in the second case B the existence of three resonances at (1.275, -0.092), (1.525, -0.038) and 
(2.011. -0.101) GeV will be used as a model constraint. We will subsequently discuss both 
cases. 

In both cases we use the formulae given in Section II. The S-submatrix for D- 
waves is then described by the nine parameters /32, /^s, Ai, A2, A3, A12, A13 and A23. The 
requirement that the S-matrix must have poles at the complex energies corresponding to 
three resonances, leads to six constraints. However, it turns out that neither case A nor 
case B allows a solution within our model if we demand that all three resonances are on the 
standard sheet for tttt, KK and pp (i.e. where all relative channel momenta have imaginary 
parts close to the physical region). On the other hand solutions for both cases are found, if 
we locate the second resonance /2(1525) on the tttt sheet, where the relative tttt momentum 
has a positive imaginary part. A physical interpretation of this would be that the coupling 
of the /2(1525) resonance to the vrvr channel is weak, meaning that /2(1525) has a rather 
small decay width into vrvr. The Particle Data Group gives that /2(1525) decays for only 
0.8% into vrvr while for 88.8% into KK, which would justify to put the /2(1525) pole not 
on the standard vrvr sheet near the physical region, but place it on the alternative vrvr sheet, 
where Im Ptttt > 0. Proceeding, the model gives results for case A as well as for case B closest 
to the data if the third resonance is put on the non-standard pp sheet. 

Then for case A, which contains the /2(1710) resonance, the parameter set is given in 



the first line of Table [V, and the model predictions for the phase shift and inelasticity are 



given by the dashed line in Fig. |^ and Fig. ^ respectively. For case A a structure near 1.6 
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GeV is present in the phase shift, and the inelasticity is much too low in comparison to the 
data. There is no evidence for a similar structure in the experimental phase shift. 

On the other hand, if we consider case B where there is no resonance at 1.7 GeV, but 
instead a resonance near 2.0 GeV, the parameter set is given on the second line of Table 
and the corresponding model prediction for phase shift and inelasticity is represented by 
the solid curves in Fig. § and Fig. |[ In the phase shift of case B the structure near 1.6 
GeV has practically disappeared, and the prediction for the phase shift (see Fig. 5) is in 
closer agreement with the data. At the same time one finds that for case B the predicted 
inelasticity (see Fig. 6) is in much better agreement with the experimental inelasticity. 

It is tempting to conclude that the presence of a /2(1710) resonance would cause more 
structure in the / = 2 phase shift than has been observed experimentally, and a value of J 
= 2 for /j(1710) is therefore unlikely. 

VI. RESULTS FOR J'^ = 1, / = 1 

Following again the formalism of Section II, the resonances imposed for / = 1 are p(770), 
p(1450) and p(1700). They occur, according to the most recent compilation of the Particle 
Data Group, at the complex energies (0.768, -0.075), (1.465, -0.,155) (1.700, -0.118) |I[ 
in units of GeV. Therefore the existence of these three J = 1 resonances again imposes 
the constraints that the analytic expression for the S-matrix for / = 1 has three poles at 
these three complex energies. These constraints determine the possible values of the nine 
model parameters, and lead to corresponding P-wave phase shifts and inelasticities. The 
parameter values for our model are given in the first line of Table |V[ The resulting phase 
shifts and inelasticities are shown respectively in Fig. ^ and Fig. |^ by the dashed curve. The 
different sets of experimental data for the / = 1 phase shift are in quite good agreement with 
each other. One observes that the model phase shifts are in good agreement with the data 
only below 1.3 GeV, but are too high above 1.3 GeV. The model phase shift also shows a 
structure near 1.6 GeV that is not present in the data. For the / = 1 inelasticity there is 
no agreement between the different experimental data sets. In one experimental set of data 
there is considerable inelasticity when the KK channel opens at 1 GeV, while in others the 
inelasticity is mainly driven by the opening of the pp channel near 1.4-1.5 GeV. Hence it is 
not possible to draw any conclusions about the model predictions from the inelasticity. The 
model result as seen in the dashed curve in Fig. ^ has only a contribution to the inelasticity 
beyond the 1.4 GeV region. It is somewhat of a surprise that a treatment of p(770), p(1450) 
and p(1700) as being very similar type resonances in all three channels, causes the model to 
fail for energies above 1.3 GeV. 

We have explored several of the options of placing one of the resonances on another sheet. 
The best agreement with the experimental data is obtained if we locate the third resonance, 
p(1700), on the non-standard KK sheet (i.e. Im pj^j^ > 0). For this case the corresponding 
parameter set is given in the second line of Table 0, and the corresponding prediction of the 
phase shift and inelasticity is given as the solid curves in Fig. |^ and Fig. ^. Comparing with 
the dashed curve, a dramatic improvement has been obtained for the phase shift. A structure 
near 1.6 - 1.7 GeV remains however present in the model. There may be an indication of 
some structure in one set of experimental phase shifts, but it is far from compelling. The 
corresponding model prediction for the inelasticity includes for this case a larger effect of 
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the KK channel, as shown by the sohd curve in Fig. ^ As long as the experimental data 
for the inelasticity suffer from the considerable internal disagreement as shown in Fig. it 
is hard to draw conclusive information from a comparison of this observable at this time. 
Other choices of possible sheets do not lead to further improvement of the predictions. 



VII. RESULTS FOR r = 2,~ , I = I 

For J'^ = 3~ there is one well established resonance p3(1690), at (1.691, -0.080) GeV, and 
a brief mentioning of p3(2250). Since we are interested in the phase shift below 2.0 GeV, 
the precise location of higher poles seems not very important. Within our model, having 
the essential ingredient of three channels and therefore nine parameters, a single resonance 
would allow too much freedom in the parameters. Just for convenience we impose therefore 
the condition that there are three resonances, even for J'^ = 3~, at energies (1.691, -0.080), 
(2.250, -0.125), (2.700, -0.300) GeV. In that case it is straightforward to obtain a good fit to 



the experimental phase shift and inelasticity by the parameter set given in Table 0. The 
model prediction for that set is given as the solid curves in Fig. |^ and Fig. 



VIII. DISCUSSION AND CONCLUSION 

In conclusion, we have constructed a three-channel model which gives a reasonable de- 
scription of the S, P, D and F-wave vrvr scattering in the 0.3 - 1.7 GeV region. Apart from 
the vrvr and KK channels, the multi-pion channels are described as effective pp and uju 
channels. As a result the model exhibits branch cuts at the KK and ppiuju) thresholds. 
Subsequently we require the S-matrix to have poles at the three lowest known resonances 
for each angular momentum J^{I) as listed by the Particle Data Group The position 
of the resonances on the various complex momentum sheets reflects the decay properties of 
each resonance. 

As a by-product of this investigation one concludes that the /o(1500) resonance plays 
a role in S-wave vr+vr^ scattering that is quite different from that of /o(980) and /o(1370). 
Also the S-wave and D-wave vr+vr^ scattering data seem to be in better agreement with a 
value of J = for /j(1710) and at the same time in disagreement with a value of J = 2. 
A further surprise for the J^{I) = 1~(1) meson resonances is that not all three resonances 
p(770), p(1450), p(1700), play a similar role in all three channels of this model. 

In order to establish how model dependent these conclusions are, it would be interesting 
to study the role of the various resonances off all angular momenta in other approaches . 
The present model can also be extremely useful to describe the final state interaction in 



reactions like ^ vr+vr [16| 
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APPENDIX A: Integrals for X'(s) 

In order to evaluate the integrals for s) of Eq. @, using Eq. ( |IOD and ([TI|) , i.e. 

X\p, s) = - dp'p' \^,^ (24) 



for real p, we define a function Z\p^ s) by 

Z\p,s) = lJdPP \p.lf^2yp._l2_,, ^ (25) 

also for real values of p. The values of Z\p,s) for / > 1 can be determined, using the 
iteration method of the Appendix of Ref. p 



Z'^\P^s) = -\^^Z\p,s). (26) 
All expressions for s) can then be determined from 

Z\p,s) = -^. (27) 
p — ip 

The expression of Eq. (^) for real p can again be analytically continued for complex p. 
The expressions for X\p, s) follow from the binomial forms 

X\p,s) = Z\p,s), (28) 

X\p,s) = Z\p,s)-(3^Z\p,s), (29) 

X\p, s) = Z\p, s) - 2f3'Z'{p, s) + (3^Z%p, s), (30) 

X'ip, s) = Z\p, s) - 3f3'Z\p, s) + 3/3^Z'ip, s) - (3'Z'{p, s). (31) 



APPENDIX B: Formulae with a second term in the tttt potential for the I = 
scattering. 

We give here explicit formulae for the Jost function and for the scattering length when 
we add, as in Ref. [2], a second term in the tttt channel for the I = scattering. The 
potential as well as the T-matrix of Eq. (P are then 4x4 matrices. We choose the index to 
label this new term in the tttt potential of range /5o and strength Aqo- We do not introduce 
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any coupling between this term and the KK and pp channels, i.e. Aoj = for z = 2, 3 and 
furthermore Aqi = 0. 

The new Jost function can then be obtained in a similar way to the derivation performed 
in section II. One gets 



D{piP2,P3) 



1 



Ao 

Rn 



Q{P2,P3) 



1 



Ao ( /3i-/3o ' 



Here Aq = Aoo/(2/3^) and Rq = l/(2/5^Xo) with, as in Eq. (|13|), 

1 



2/3oi/3o-^Ply' 



Eq. (^2|) contains the expressions 



n( ^ 1 ^ ^2 ^ ^3 ^ A2A3 - Ajg 

Q(P2,P3) = ^ + -B- + -B- + -B-B : 

K2 Ks -K2-K3 



P{P2,P3) 

Ri 



(32) 



(33) 



(34) 



P(P2,P3) = 1 + 



A1A2-A22 AiAg-Afg 



R2 



R3 



-^123 
R2R3 



with 



(35) 



^123 



2A12A23A13 - A^As - A?2A3 - A^gAi + A1A2A3. 



It can be seen that if Aq = 0, Eq. ( |32D reduces to Eq. ( pT]) as it should. 
The TTTT scattering length can be calculated as 



lim 



As limpj^o^i = 1; one obtains 



lim 



D{-P1,P2,P3) - D{pi,P2,P3) 

2ipiD{pi,p2,P3) 



(36) 



(37) 



-2P{p2,P3)/Pi-2Ao/Po 


Q{P2,P3)+P{P2,P3)^1^\ 


<5(P2,P3)(l + Ao)+P(p2,P3) 




2 





(38) 



In Eq. (|38| ) = ^1 ~ ""^2 ^'^'^ P3 = "^1 ~ 
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FIGURES 




FIG. 1. S-wave tttt phase shifts from the three-channel model, if poles are present for /o(980), 
/o(1370) and /o(1500). The dashed curve is the result for Z = if the /o(1500) pole occurs at 
Im Ptttt < with the parameter set of the first line in Table |lll The solid curve is the prediction 



for the parameter set in the second line in Table III, where the /o(1500) resonance pole corresponds 
to Im Ptttt > 0. Data are from Refs. [pT,l^,l^. 
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FIG. 2. S-wave vrvr inelasticities from the three-channel model, if poles are present for /o(980), 
/o(1370) and /o(1500). The dashed curve is the result for / = if /o(1500) has Im Ptttt < and 
the parameters of the first line in Table p^ . Solid curve is the prediction for the parameter set in 
the second line in Table III, where /o(1500) has Im Ptttt > 0. Data are from Refs. [11,12,14]. 
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FIG. 3. S-wave vrvr phase shifts from the three-channel model, if poles are present for /o(980), 
/o(1370) and /o(1710). Solid, dashed and dash-dotted curves are the prediction for / = parameter 
sets from Table lines three thru five. Data are as in Fig. |l[ 
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3 CHANNELS; L = 



<; 



1 ' ' ' 1 1 — 

S-WAVE 77-77 SCATTERING INELASTICITY 

INPUT 3 POLES: fo(980), f(,(1370), fo(1710) 




M(pi-pi) (GeV) 



FIG. 4. S-wave vrvr inelasticities from the three-channel model, if poles are present for /o(980), 
/o(1370) and /o(1710). Curves are as in Fig. ^. Data are as in Fig. ^. 
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FIG. 5. D-wave vrvr phase shifts from the three-channel model, if poles are present for /2(1270), 
/2(1525) and /2(2010). Dashed curve is the prediction from the first / = 2 parameter set from 
Table |l^, case A, with the /2(1525) resonance with Im p^^.,^ > 0, while the third resonance is the 
/2(1710). Solid curve is the prediction from the second / = 2 parameter set from Table |lV| case 
B, with the /2(1525) resonance with Im p^^-^ > 0, while the third resonance is the /2(2010). Data 
are from Ref. []l2| . 
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3 CHANNELS; L = 2 pi-pi INELASTICITY 



1 1 1 1 1 1 1 1 1 1 1- 

SOUD CURVE FOR INPUT: tg(1270), fg(1525), fg(2010) 
DASHED CURVE IF INPUT: f2(1270), f3(1525), fg(1710) 



<; 
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M(pi-pi) (GeV) 



FIG. 6. D-wave vrvr inelasticities from the three-channel model, if poles are present for /2(1270), 
/2(1525) and /2(2010). Dashed and solid curves are as in Fig. ^. Data are from Ref. p^ . 
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3 CHANNELS; L 
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FIG. 7. P-wave vrvr phase shifts (in degrees) from various experimental analyses. The curves 
represent vrvr phase shifts for / = 1 from the three-channel model, if poles are present for yo(770), 
/9(1450) and p(1700). The dashed curve has all resonances on the standard vrvr, KK and pp sheets, 
while the solid curve is the result if p(1700) is located on the sheet where Im p^k > 0- Data points 
are from Refs. |1C,12|. 
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FIG. 8. P-wave tttt inelasticities from various experimental analyses. The curves represent tttt 
inelasticities for / = 1 from the three-channel model, if poles are present for p{770), p(1450) and 
/9(1700). The dashed curve has all resonances on the standard vrvr, KK and pp sheets, while the 
solid curve is the result if yo(1700) is located on the sheet where Im p^k > 0- Data points are from 
Ref. |10|,||. 
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3 CHANNELS; L = 3 
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FIG. 9. F-wave vrvr phase shifts from the three-channel model, if poles are present for /93(1690), 
/93(2250) and p3(2700). Solid curve is the prediction for the / = 3 parameter set from Table [V^ . 
Data are from Refs. p^ , [T3| . 
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3 CHANNELS; L 



1 1 1 1 1 1 1 1 1 — 

F-WAVE 77-77 SCATTERING INELASTICITY 
INPUT 3 POLES: P3(1690), Pg(Z25Q), /Og(2700) 
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FIG. 10. F-wave vrvr inelasticities from the three-channel model, if poles are present for 
yO3(1690), /93(2250) and p3(2700). Solid curve is the prediction for the / = 3 parameter set from 
Table 0. Data are from Refs. HJll. 
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TABLES 

TABLE L Relations of Ajj to Xij and Ri to Xi for angular momentum /. 



l_ ^ Ri 

A,,/(2/33) A2./(4/33/33) l/{2(5lX,) 

1 A,,/(16/?f) A2./(256/3f/3|) l/(16/3fX,) 

2 A,./(256/3j) A2,./(2562/3//3j) l/(256/3jX0 

3 A2./(20482/39/?9) l/(2048/3f X,) 



TABLE IL Panicle Dala resonances considered in lliis work. 



r{I) 0+(0) 2+(0) 3-(l) 

~^ /o(980) p(770) hivm) P3(1690) 

— /o(1370) p(1450) /2(1525) P3(2250) 

— /o(1500) /9(1700) /2(1710)? /93(2700) 

— /o(1710)? p(2150) /2(2010) - 



TABLE III. 


DifCereut 


model 


]")aram('t e 


'rs for / 


= 0. i.e. 


ni) = 0^ 


-(0) (S in 


Gi-V). Figs 


. 1 1. 




Pi 


(h 




Ao 


Ai 


As 


A3 


A12 


Al3 


A23 


2.80 


0.85 


0.40 


1.75 


-2.30 


4.05 


-0.158 


-1.300 


0.5523 


0.0399 


0.5425 


3.60 


1.35 


0.50 


0.25 


-2.70 


3.95 


-0.678 


-5.561 


0.3059 


0.3173 


2.0061 


3.80 


1.25 


0.20 


1.00 


-2.40 


3.80 


-0.946 


-1.575 


0.1689 


0.3487 


1.2793 


2.80 


0.85 


0.30 


3.25 


-2.30 


4.20 


-0.124 


-1.013 


0.5646 


0.3230 


0.4532 


4.00 


1.25 


0.30 


0.50 


-2.30 


3.70 


-0.914 


-2.656 


0.1333 


0.5098 


1.6397 



TABLE IV. Different model parameters for / = 2, i.e. J'^{I) = 2+(0) (/3 in GeV), Figs. 5-6. 



Pi 


f32 




Ai 


A2 


A3 


A12 


Al3 


A23 


2.400 


1.900 


0.750 


-0.210 


-0.252 


-0.050 


0.0428 


0.1222 


0.0095 


3.600 


0.300 


1.750 


-0.300 


1.596 


-0.321 


0.1059 


0.0107 


0.0977 



TABLE V. Different model parameters for / = 1, i.e. J''(/) = 1-(1) (/3 in GeV), Figs. 7-8. 



Pi 


P2 


P3 


Ai 


A2 


A3 


A12 


Ai3 


A23 


3.200 


15.500 


3.250 


-0.855 


-0.991 


-1.051 


0.0035 


0.0843 


0.0074 


2.600 


0.500 


0.750 


-0.560 


2.966 


-1.457 


0.9756 


0.4386 


0.2775 
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TABLE VI. Model parameters for / = 3, i.e. = 3-(l) (/? in GeV), Figs. 9 - 10. 



Pi 


/?2 




Ai 


A2 


A3 


A12 


Al3 


A23 


3.800 


6.500 


0.750 


-0.170 


-0.188 


0.048 


0.0013 


0.0413 


0.0146 
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